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Abstract 

We study the Euler-Poisson equations of describing the evolution of the 
gaseous star in astrophysics. Firstly, we construct a family of analytical 
blowup solutions for the isothermal case in R^. Furthermore the blowup 
rate of the above solutions is also studied and some remarks about the ap- 
plicability of such solutions to the Navier-Stokes-Poisson equations and the 
drift-diffusion model in semiconductors are included. Finally, for the isother- 
mal case (7 = 1), the result of Makino and Perthame for the tame solutions 
is extended to show that the life span of such solutions must be finite if the 
initial data is with compact support. 

1 Introduction 

The evolution of a self-gravitating fluid such as gaseous stars can be formulated by 
the Euler-Poisson equation of the following form: 

Pt+V ■ (pu) =0, 
(pn)^+V ■ (pw ® u) + VP =-pV<l', , . 

St + u-VS = 0, ^ ' 

A^t,x) = a{N)gp, 

where a{N) is a constant related to the unit ball in R^: a(l) = 2; a{2) = 2tt; For 
simplicity, we take the constant term g = 1. For > 3, 

a{N) = N{N - 2)V{N) = N{N - 2) ^^^^7+1) ' 
*E-mail address: nevetsyuen@hotmail.com 
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where V^(A^) is the volume of the unit ball in i?^ and F is a Gamma function. As 
usual, p = p(t,x), u = u(t,x) G and S(t,x) are the density, the velocity and 
the entropy respectively. P = P{p) is the pressure. In the above system, the self- 
gravitational potential field $ = x) is determined by the density p through the 
Poisson equation. The equations ([T])i and ([TJ2 are the compressible Euler equation 
with forcing term. The equation ([T])3 is the Poisson equation through which the 
gravitational potential is determined by the density distribution of the gas itself. 
Thus, we called the system ([T]) the Euler-Poisson equations. Here, the viscosity term 
does not appear, that is, the viscous effect is neglected. In this case, the equations 
can be viewed as a prefect gas model. For = 3, ([1]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [2], P| for a detail about the system. 

If we take S(t, x) = In K, for some fixed K > 0, we have a 7-law on the pressure 
Pip), i.e. 

P{p) = Kp^ = l (2) 

7 

which is a commonly the hypothesis. The constant 7 = Cp/c^, > 1, where Cp, Cy 
are the specific heats per unit mass under constant pressure and constant volume 
respectively, is the ratio of the specific heats, that is, the adiabatic exponent in ([2]). 
In particular, the fluid is called isothermal if 7 = 1. It can be used for constructing 
models with non-degenerate isothermal cores, which have a role in connection with 
the so-called Schonberg-Chandrasekhar limit [S]. 

For the physical dimension = 3, we are interested in the hydrostatic equilib- 
rium specified by u = 0, 5* = In A'. According to [2], the ratio between the core 
density p(0) and the mean density p for 6/5 < 7 < 2 is given by 

-y {z) 



z=zo 



p(o) V z 

where y is the solution of the Lane-Emden equation with n = 1/(7 — 1), 
y+-y + y^ = 0, y{0) = a > 0, y(0) = 0, n ^ 



z 7 — 1 

and zq is the first zero of y{zo) = 0. We can solve the Lane-Emden equation 
analytically for 



yanal \'^J 



l-\z\ n = 0; 



, \ , n = 5, 



and for the other values, only numerical values can be obtained. It can be shown 
that for n < 5, the radius of polytropic models is finite; for n > 5, the radius is 
infinite. 
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Gambin [6J and Bezard [Ij obtained the existence results about the exphcitly 
stationary solution (m = 0) for 7 = 6/5 : 

,.f?^V^,..V)-. (3) 



V 27r 

The Poisson equation can be solved as 



^{t,x)= G{x -y)p{t,y)dy, 

where G is the Green's function for the Poisson equation in the A^-dimensional spaces 
defined by 

f kl, N = l- 

G{x) = I log\x\,N = 2; 

iV>3. 

In the following, we always seek solutions in spherical symmetry. Thus, the Poisson 
equation ([TJa is transformed to 

r^-'^rr (t, x) + {N- 1) r^-2$,,=a (A^) pr^"\ 
^ Jo 

Definition 1 (Blowup) We say a solution blows up if one of the following condi- 
tions is satisfied: 

(1) The solution becomes infinitely large at some point x and some finite time T; 

(2) The derivative of the solution becomes infinitely large at some point x and some 
finite time T. 

In this paper, we concern about blowup solutions for the A^- dimensional isother- 
mal Euler-Poisson equations, which may describe the phenomenon called the core 
collapsing in physics. And our aim is to construct a family of such blowup solutions 
to it. 

Historically in astrophysics, Goldreich and Weber ^ constructed the analytical 
blowup solution (collapsing) solution of the 3-dimensional Euler-Poisson equation 
for 7 = 4/3 for the non- rotating gas spheres. After that, Makino ^ obtained the 
rigorously mathematical proof of the existence of such kind of blowup solutions. 
And in [5j, we find the extension of the above blowup solutions to the case A^ > 3 
and 7 = {2N — 2)/N. In [11], the solutions with a from is written as 



p{t,r) = { ^itry^aUf^^''''''^ for r < ait)Z„ ^ 
' ^ ' 0, for a{t)Z^, <r. ' V , ; a{t) ' v , ; 

~ a(t)^-i ' ^(0) = «o > 0' «(0)= Oi 



y{z)+^y{z)+^^^^=^y{zf/^^-'^= p, y{0) = a > 0, y{0) = 0, 

(4) 
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where /z = [N{N - 2)\]/{2N - 2)K and the finite is the first zero of y{z). 
And no other analytical blowup solution for the Euler-Poisson equation has been 
obtained. Through there are a lot of numerical simulation results on the Euler- 
Poisson equation. 

In section 2, we obtain the blowup solutions for the Euler-Poisson equation in 
spherical symmetry in the 2-dimensional case, 



in the form of the following theorem. 

Theorem 2 For the 2-dimensional Euler-Poisson equations in radial symmetry, 
there exists a family of solutions. 



where K > 0, p = 2X/K with a sufficiently small A and a are constants. 

(1) When A > 0, the solutions blow up in a finite time T; 

(2) When A = 0, if ai < 0, the solutions blow up at t = — ao/^i- 

The blowup rate of the solution is discussed in section 3. Actually, we have 
Theorem 3 The blowup rate of the solution ^ is. 



Pt + upr + pUr + }.pu = 

p {Ut + UUr) + Kpr = - 



f^J^p{t,s)sds, 



(5) 



p{t,r) = 
a{t)= - 

yix)+^yix)^ 




(6) 



limp(t,0) {T-tr>0{l), 



with rj < 2. 



In the last section, we can rewrite ([T]) in scalar form. 



9p I spN dp_ spN duk 



1,2, ...AT. 



(7) 



The tame solution is introduced in here. 



Definition 4 A solution {p,u) for ^ is called "tame" solution if 
(l){p,u) G C^{[0,T) X R^), p > 0, p{t) has compact support, and 
^^;p(7-i)/2 ^ C\[0,T) xR^),ue C([0,T);50), and 



dui sr^N dui (9$ , 



8 



holds on the exterior of the support of p. 
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If only the condition (1) is assumed, the solution is called "classical". For 7 > 1, 
under the transformation 

w = —p^ ' " , 

7 - 1 

Makino and Perthame [10] showed that the life span of tame solution of the Euler- 
Poisson equation in spherical symmetry is finite. As we are interested in non-global 
existence for the isothermal case (7 = 1) in 2 dimension, i.e. 

pt + upr + pUr + \pu = 0, 

p {Ut + UUr) + Kpr = /J' p(t, S)sds, 

the result of Makino and Perthame [10] for the tame solutions is extended as the 
following theorem. 

Theorem 5 Let {p{t),u{t)) be a radially symmetric tame solution of ^ on < 
t < T. If the support of (p(0),u(0)) is compact and p(0) is not identically equal to 
zero, then T must be finite. 



2 Separable Blowup Solutions 

In this section, before presenting the proof of Theorem [2], we prepare some lemmas. 
First, we obtain a general class of solutions for the continuity equation of mass in 
radial symmetry (I5])i. 

Lemma 6 For the equation of conservation of mass in radial symmetry 



there exist solutions, 



1 , , 

Pt + UPr + pUr + -pU = 0, (10) 

r 



fir /ait)) , , ait) 
a(t)^ ait) 



with the form with f > E and a{t) > G 

Proof. We just plug ^ into ([10]). Then 

1 

Pt + UPr + pUr + -pU 

r 

-2a{t)f{r/a{t)) a{t)rf{r/a{t)) 



a(t)3 a(t)4 

a{t)rfir/ait)) f jr / a{t)) ajt) 1 f {r / a{t)) a{t) ^ 

a{t) a(t)3 a(t)2 a{t) r a{t)^ a{t)' 
= 0. 



6 



M.W.Yuen 



The proof is completed. ■ 

Secondly, we obtain an estimate to the equation, 

«W = -4)' 

a(0) = ao > 0, a(0) = ai. 



(12) 



Lemma 7 For the Emden equation u^) . we have, 

(1) If X > 0, there exists a finite time T_ < +oo such that a{TJ) = 0; 

(2) If X = 0, it holds that for any t > 0, 

a{t) = ao + ait. 
Proof. (1) By integrating ([T2l) . we have 

<-a{tf = -X\na{t) + e (13) 
2 

where 9 = Xlnao + ^a\. 
From f|T3l) . we get, 

ait) < e'/\ 
If the statement (1) is not true, we have 

< a{t) < e^/^, for all t > 0. 

But since 

A -A 
"^'^ = - ^' 

we integrate this twice to deduce 

J ^dsdr + Cit + Co = + Cit + Co- 

By taking t large enough, we get 

a{t) < 0. 

As a contradiction is met, the statement (1) is true. It is trivial to verify (2). ■ 

Remark 8 The modified systems of the above Emden equation may he referred more 
detailedly m Jf, W^,^^. 

Lemma 9 There exists a sufficiently small xq > 0, such that the equation 

y{x)+^y{x)+(rey^^'^= fi, 
1/(0) = «, y{0) = 0, 

where a > 0, ^, and a are constants, has a solution y = y{x,^) G C^[0,xo]. 
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Proof. The lemma can be proved by the fixed point theorem. Multiply (fT4|) by x, 
to give 



Notice y{0) = 0, we have 



y{x) 



By using y{0) = a, IHM is reduced to 



y{x) 



X 



- aey^'^)ds, y{0) = a. 



Set 



f{x,y{x)) =- [ s{ii-aey^''^)ds. 



X 



then for any Xq > 0, we get / G C^[0, Xq]. and for any yi^ y2 G C^[0, Xq], we have, 



\f{x,yi{x)) - f{x,y2{x) 



X 



As is a function of y, we can show that the function e^, is Lipschitz-continuous. 
And we get, 

\f{x,yi{x)) - f{x,y2{x)\ 
0{l)J;;s\{y2{s)-y,is)\ds 

X 

< 0{l)xo sup \yi{s) -i/2(s)| , 

0<x<xo 

where r G [0, x] C [0, Xq]. Let 

/■X 

Ty{x) = a+ f{s,y{s))ds. 
Jo 

We have Ty G C[0, xq] and 

\Tyi{x) -Ty2{x)\ 

f{s,yi{s))ds - / f{s,y2{s))ds 



< O(l)xo sup \y{x)i - y{x)2\ . 

0<x<xi) 

By choosing xq > to be a sufficiently small number, such that O(l)xo < 1, this 
shows that the mapping T : C[0, Xq] — > C[0, Xq], is a contraction with the sup- 
norm. By the fixed point theorem, there exists a unique y{x) G C[0, Xq], such that 
Ty{x) = y{x). The proof is completed. ■ 
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Lemma 10 The equation, 

y{x) + \y{x)+aey^-^= 0, 
y(0) = a, y(0) = 0, 

where a > and a are constants, has a solution in [0, +oo) and lim y{x) = — oo. 
Proof. By integrating ( fTSl) . we have, 



a;— >+oo 



X 



y(x) = -- I sey'^'^ds < 0. (16) 
X Jo 

Thus, for < X < Xq, y{x) has a uniform lower upper bound 

y{x) < y{0) = a. 

As we obtained he local existence in Lemma [H there are two possibilities: 

{l)y{x) only exists in some finite interval [0,Xo]: (la) lim y{x) = — oo; [lh)y{x) has 

an uniformly lower bound, i.e. y{x) > ao for some constant ao- 

{2)y{x) exists in [0, +oo): (2a) lim y{x) = —oo; {2h)y{x) has an uniformly lower 

bound, i.e. y{x) > a for some constant oq. 

We claim that possibility (1) doesn't exist. We need to reject (lb) first: If the 
statement (lb) is true, f|T6|) becomes 

rx 

- axe"" = / xe'^ds < y{x). (17) 



X 







Thus, y{x) is bounded in [0,Xo]. Therefore, we can use the fixed point theorem 
again to obtain a large domain of existence, such that [0, Xq + 6] for some positive 
number S. There is a contradiction. Therefore, (lb) is rejected. 
Next, we do not accept (la) because of the following reason: It is impossible that 
lim y{x) = —oo, as from f|T7|) . y{x) has a lower bound in [0, Xq]: 



crxoe" < y{x). 



Thus, ffTSl) becomes. 



PXO 

yixo) = yiO) + / yix)dx 
Jo 

> a — axoe^dx 
Jo 

= a — axle". 
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Since y{x) is bounded below in [0, xq], it contracts the statement (la), such that 
lim y{x) = — oo. So, we can exclude the possibility (1). 

We claim that the possibility (2b) doesn't exist. It is because 



y{x) 



a 

-— I se 

X 



yi')ds <-- I e^'^sds 



X 



Then, we have. 



y{x) <a- 



ae 



ao 



-X 



(19) 



By letting x — *■ oo, f[T9|) turns out to be, 

y{x) = -oo. 

Since a contradiction is established, we exclude the possibility (2b). Thus, the 
equation f|T5|) exists in [0, +oo) and lim y{x) = — oo. This completes the proof. ■ 

x—*+oo 

Now transferring f|T^ to the first-order system, 

dy 



dy 
dx 



dx 

1, 



- aey + /i, 

we consider the system along the solution curve y = y{x,0), y = y{x,0), < xq < 
X < +00. Since the right-hand side of this system is continuously differentiable in 
X, y, y and /i. We apply the comparison theorem (Theorem 7.4) in [3], to obtain 
the following lemma. 

Lemma 11 For a sufficiently small e such that \fi\ < e, the equation, 

y{x)+ly{x)+crey^''^= /i, 
y(0) = a,yiO)=0, 

exists in [0, oo), and lim y{x) = — oo. 



X— > + oo 



Now, we give the proof of Theorem [21 
Proof of Theorem 2. From Lemma [TTl we easily get that satisfy ([HDi. For the 
momentum equation ([5]) 2, we get. 



r 

p{Ut + UUr) + Kpr H I p{t, s)sds 

r J 



P 



a{t) 

P 
a{t) 

P 



a{t)r+Ky{ 



2ti 



a{t) ra{t) 



M^l<t))sds 



Xr ■ , r , 2n 

+ Ky{^) + 



r/a{t) 



a{t) 
r 



'a{t)' 



■■ait)' 
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Here, we use the property of a{t): 

A 



a{t)-- 



a(t)' 
and denote 

X 

Qi^) = Qi^) = -^x+Ky{x)+— [ ey'-'^^sds. 
a{t) X J 



Differentiate Q{x) with respect to x, 

X 

Q(x) = -A+ir?/(x)+27re^(^) / e^^'^sc/s 

X J 


X 





X 



j Ax + Ky{x) - Kfix + j e^^'hds 



= --Q(x). 

X 

The above result is due to the fact that we choose the following ordinary differential 
equation: 

1/(0) = a, y(0)=0. 

With Q{0) = 0, this implies that Q{x) = 0. By using Lemma [7] about a{t) and 
Lemma [10] about y{x), we have shown that the family of the solutions blows up 
in finite time T under the prescribed conditions of Theorem [21 This completes the 
proof. ■ 

Remark 12 Besides, the above solution requiring y{0) > is not necessary for ^ 
in contrast to in f^, f^. And the mass of the above solution is 



r,„ r+OO 



M=l p(t,s)ds= I 4^e^^'^"^°))srfs = 27r / e^^'^sds, 
Jr^ Jo ai^r Jo 

where ri{N) is the measure of a unit ball. Thus the mass of the solution ^ depends 
on the initial data 1/(0). And it is not easy to determine if the mass is finite or 
infinite with different datum. It is different for the finite mass of the family of 
solutions for 7 = 4/3 and N = 3 which is independent of y{0) in J^. 

Remark 13 Now, we consider the stationary solutions of (J^), in spherical symme- 
try, i.e. 

Pr = -p<l>r. (20) 
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For 7 = 1, nW^I becomes, 

p = e-*/^+^, 

where C is a constant. 

For our convenience, C = is chosen, 

p = e-*/^. 

From the Poisson equation (QPs? obtain, 

A$ = a(Ar)e-*/^. (21) 
As the solution is in radially symmetric, [21\) becomes, 

<^rr + -$r - 27re"*/^ = 0. 



■ r 

r 



Under the transformation. 



y{x) = 



we have 



K ' 

1 27r y''^^ 

y[x)+-y{x) +—e = 0. (22) 
X K 

Furthermore by setting A = and a(0) = in Theorem [21 we easily have the 
following corollary for the non-trivial stationary solutions. 

Corollary 14 For the 2-dimensional isothermal Euler-Poisson equation, there ex- 
ists a family of stationary solution, i.e. 

p{t,r) = ^e ■ 
y(x)+iy(x)+f e^^")= p, yiO) = a, y(0) = 0, 

where a is an arbitrary positive constant. 

Remark 15 Our blowup solutions only work for the 2-dimensional case. But we do 
not know what will happen after the critical time that the solutions blow up. 
If we consider the system of conservation laws with viscosity, i.e. with u > 0, 

Pt+V ■ (pu) =0, 
{pu)^+V ■ {pu (g) m) + VP =-pV$ + uAu, 
A^{t,x) = 27rp. 

We get the corresponding Navier-Stokes-Poisson equation. The equation describes 
the situation closer to the real model compared with the Euler-Poisson equations in 
gaseous stars. And our family of solutions is also suitable for it. This is due to the 
fact in spherical symmetry for the vector Laplacian in u{t,r): 

1 1 

lAU = Urr H Mr ^W. 
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Remark 16 In /H] /. we have the lemma to control the modified Emden equation 
with (3 > 0, 

ait) + Pa(t)= — , a(0) = a„ > 0, a(0)= ai. 

a{t) 

The similar blowup results may be obtained for the 2-dimensional Euler-Poisson 
equation with frictional damping or Navier-Stokes-Poisson equation with frictional 
damping or not, i.e. 

Pt+V ■ (pu) =0, 
(pn)j+V • {pu (g) u) + VP = -pV$ - Ppu + uAu, 
A^{t,x) = 2np, 

where /3 > and u >0. The corresponding family of blowup solution is 

p(t,r) = 4,e^('-H*)),«(t,r) =|§r; 
a{t) + Pa{t)= - a(0) = > 0,a(0)= oi; (24) 
y{x)+ly{x) + ^^e'^^^= p,y{0) = a,y{0) = 0, 

where p = 2X/K and a are constants. 

Remark 17 Besides, if we consider the drift- diffusion model in semiconductors, 

Pt+W ■ (pu) =0, 
(pu)i+V • {pu O m) + VP = +pV$ - ppu + uAu, 
A$(t,x) = 27rp, 

the special solutions with infinite mass may be obtained as follows 

p{t,r) = ^ey^^/-('^\u{t,r) =|gr; 
a{t) + I3a{t)= - a(0) = %> 0,a(0)= oi; (25) 



2{t) ^ 



yix)+^yix) - f e"^"^= p,i/(0) = a,j/(0) = 0. 



3 Blowup Rate 

In this section, we present the confirmation of Theorem [31 The blowup rate of the 
constructed solutions (EI), for the Euler-Poisson equation is studied. It is interesting 
to investigate how fast the blowup solution tends to infinity as the time tends to the 
critical value T. 

Proof of Theorem [Si We choose a finite time t to make a{t) sufficiently small 
enough, such that —A lna + 6 > 0. From the lemma[71 it is clear for that there exists 
a finite time to > such that a(t) < for t > to. Next, since 
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when t > to + e = ti, this means there exists a constant e > 0, such that 

a{t) < e'/^ - 5, 

where 5 is a sufficiently small positive constant, we have a smaller upper bound of 
a{t). On the other hand, as 



a = -V-2Alna + 6. 

It is easy to see that, 

f' , r^*^ da 0(l)da 

t = ti+ ds = ti- ^ >ti+ / ' . (26) 

Jh V-2A Ina + ^ Jo V-\na 

We denote 5 as a sufficiently small positive number, such that x G (0, 6], 

— Inx < — ^, 

where the constant S > 0. Then, fl2B]) becomes, 

r(*) 0(l)da ^, , , ,./2+i 
Jo \/^laP 

By letting t\ — > T_, we get. 



hm ait) {T-t,)-'/^'^'^<0{l). 

As y{0) = a, we can estimate the blowup rate at the origin by 

^hmp(t,0) (T-t)^/("+^)>0(l). 



Thus, we obtain. 



^hmp(t,0) (T-t)" >0(1), 



with r] < 2, and complete the proof. 



4 Non-global Existence for Tame Solutions 

In this section, we present the proof for Theorem [51 The technique to Makino and 
Perthame's [10] is similar. 

Proof of Theorem [5] . We denote the radius of the compact support of the tame 
solution by R{t). Under the assumption that the initial data has a compact support, 
i.e. p(0, r) = u(0, r) = for r > -R(O), from (^2, we must have 



Kpr = 0. 
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And the equation ([9]) 2 becomes 



Ut + UUr = / p{t, s)sds, 



~dt' 



along the curve r = r{t; 0, y) where r = r(t; to, Tq) is the solution of the charac- 
teristic equation outside the support of the tame solution, i.e. 

— = u{t,r), r \t=to = tq. 

Clearly we have for y > R{t), 

-u(t, r(t; 0, y)) = Ut + uUj. = — $r = G(r) = / p(t, s)sds < 0. 

' r Jo 

Therefore, we get 

-rit;0,y) = uit,rit;0,y)) = G(r, r(r; 0, 2/))rfr < 0, 

as ^(O,^) = 0. Hence r{t;Q,y) < ?/ for < t < T, R{t) < y. For p{t) = outside 
the compact support of the tame solution, this implies 

G(t,r(t,0,i?(0))): ^ 



r(t,0,i?(0))' 
where 



M = I p{t,x)dx = 2TT / p{t,s)sds, 

JR^ Jo 

is the total mass and independent of t. And If p(0) is not identically equal to zero, 
as the solution belongs C^, we have M > 0. Therefore r = r{t; 0, -R(O)) satisfies the 
equation of free fall 

d'^r M 

— = --,r(0) = i?(0),r(0) = 0. (27) 

From this it follows 
This implies 



The proof is completed. ■ 
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